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Executive  Summary 


This  documentation  is  the  final  technical  report  for  work  performed  under  AFOSR  grant 
F49620-96-1-0319.  The  title  of  this  project  was  "Atomistic  and  Ab  Initio  Calculations  of  Ternary 
II-IV-V2  semiconductors"  and  the  Principal  Investigator  was  Ravi  Pandey,  Physics  Department, 
Michigan  Technological  University,  Houghton,  MI  49931.  It  was  36-months  research  program  to 
use  atomistic  and  first-principles  methods  to  study  chalcopyrites  and  related  semiconducting 
materials.  During  the  grant's  active  period,  the  following  repons  have  been  published  : 

(i)  M.  Ohmer  and  R.  Pandey  (Guest  Editors),  MRS  Bulletin  July  1998 

"Emergence  of  chalcopyrites  as  nonlinear  optical  materials. " 

(ii)  P.  Zapol,  R.  Pandey,  M.  Seel,  J.  M.  Recio  and  M.  Ohmer  and  J.  M.  Recio  1999, 

J.  Phys  :  Condensed  Matter  11,1. 

"  Density  functional  study  of  the  structure,  thermodynamics  and 
electronic  properties  of  CdGeAs2." 

(iii)  R.  Pandey.  M.  Ohmer  and  J.  Gale.  1998,  J.  Phys  :  Condensed  Matter  10,  5525. 

"A  theoretical  study  of  native  acceptors  in  CdGeAs2. " 

(iv)  R.  Pandey,  M.  Ohmer,  A.  Costales  and  J.  M.  Recio,  1998,  MRS  Symposium 

Proceedings  Volume  484,  p.525. 

"  Atomistic  calculations  of  dopant  binding  energies  in  ZnGeP2.  ” 

(v)  R.  Pandey.  M.  Ohmer.  A.  Co.stales  and  J.  M.  Recio,  1999,  MRS  Symposium 

Proceedings,  Submitted. 

"  Modeling  of  dopant  properties  in  CdGeAs2- " 

(vi)  P.  Zapol,  R.  Pandey  and  J.  Gale,  1997,  J.  Phys.  Condensed  Matter  9,  9517. 

"  An  interatomic  potential  study  of  the  properties  of  GaN. " 

(vii)  R.  Pandey,  P.  Zapol  and  M.  Cau.sa,  1997,  Phys.  Rev.  B  55,  R 16009. 

"  /i  theoretical  .study  of  nonpolar  surfaces  of  AIN. " 

(viii)  R.  Pandey.  M.  Rcrat  and  M.  Causa,  1999,  Appl.  Phys  Letts.  75,  In  Pre.ss. 

"  First-principles  study  of  stability,  hand  structure,  and  optical 
properties  of  the  ordered  Geo.soSnoso  alloy." 
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Accomplishments/New  Findings : 

(1)  Atomistic  defect  calculations  in  ZnGePi  have  identified  the  dominant  native  acceptor  to 
be  the  zinc  acceptor  vacancy  which  is  responsible  for  the  presence  of  absorption  band  near  the 
pump  wavelength  limits.  Although  cation  lattice  disorder  is  predicted  to  be  the  dominant  native 
defect,  the  dominant  defects  controlling  the  properties  of  interest  would  seem  to  be  the  zinc 
acceptor  vacancy  (binding  energy  of  0.57  eV)  which  is  partially  compensated  by  a  phosphorus 
donor  vacancy  (binding  energy  of  0.68  eV).  For  the  EPR-active  acceptor  center,  the  calculated 
binding  energies  and  lattice  distortion  corroborates  the  ENDOR  spectrum  a.s.sociating  the  zinc 
vacancy  with  the  acceptor  center. 

(2)  In  CdGeA.S2.  the  acceptor  center  identified  by  Hall-effect  measurements  and  EPR  is 
found  to  be  related  to  the  delocalized  hole  shared  by  the  four  /\.v  neighbors  bound  to  Cdpe-  For 
this  center,  calculations  yield  a  binding  energy  of  0.13  eV  in  an  agreement  with  the  experimental 
value  of  0. 15  eV  obtained  by  the  Hall-effect  measurements. 

Altlioui’li  ZnCeP2  and  CdGeAs2  heloni’  to  the  same  chalcopyrite 
family,  the  nature  of  dominant  acceptors  in  these  materials  is  predicted  to  be  different.  Defect 
calculations  corrohoratin{>  the  Hall  effect  (Ohmer)  and  magnetic  resonance  (Halliburton) 
studies  find  that  the  zinc  vacancy,  not  the  zinc  antisite  (Znce)-  as.wciated  with  the  dominant 
acceptor  center  in  ZnGeP2-  This  is  not  the  ca.se  in  CdGeAs2  where  a  cadmium  antisite  ( Cdce)  ts 
predicted  to  be  as.sociated  with  the  dominant  acceptor  in  the  lattice.  We  believe  that  this 
difference  may  well  be  due  to  defect-induced  lattice  distortion  which  plays  a  key  role  in 
stabilizing  the  hole  slates  in  the  lattice. 

(3)  ZnGeP2  is  the  NLO  material  for  mid  IR  lasers  allowing  high-power  tunability  in  the 
spectral  region  of  2-5  pm.  An  absorption  band  around  1-2  pm  is,  however,  found  to  affect  the 
usefulness  of  ZnGePz-  Experimental  efforts  involving  selective  doping  of  ZnGeP2  are  underway 
at  Lockheed  to  reduce  the  concentration  of  the  acceptor-complex  giving  rise  to  an  absorption 

band  around  1-2  pm  in  the  lattice.  To  provide  a  guidance  to  experimental  efforts,  we  have 
performed  atomistic  calculations  on  ZnGeP2  involving  cation  dopants,  namely  Cu,  Ag,  B,  Al, 
Ga  and  In.  Calculations  predict  small  binding  energies  for  Cu  and  Ag  substituting  Zn  in  the 
lattice  which  are  in  agreement  with  the  available  experimental  data.  The  group  HI  dopants  (i.e.  B, 
Al,  Ga  and  In)  at  the  Ge  site  are  predicted  to  have  large  binding  energies  for  a  hole  except  B 
which  shows  a  distinct  behavior.  This  is  due  to  large  mismatch  in  atomic  sizes  of  B  and  Ge.  At 


(iii) 


the  Zn  site,  the  calculated  binding  energies  of  the  group  III  dopants  place  donor  levels  in  the 
middle  of  the  band  gap. 

(4)  Group-IV  semiconductor  alloys  have  immen.se  potential  for  applications  in  the  next 
generation  of  Si-based  electronic  and  photonic  devices.  The  MBE  group  at  Michigan  Tech 
funded  by  DARPA  has  recently  reported  (Phys  Rev  Lett  80,  1022  (1998))  the  morphological 
evolution  of  an  epitaxially  strained  Ge-Sn  alloy  thin  films.  Since  a  knowledge  of  structural, 
electronic  and  optical  properties  of  GeSn  alloys  is  critical  in  promoting  alloy  thin  films  for 
device  applications,  we  have  performed  a  first-principles  study  on  the  ordered  Ge0.50Sn0.50  alloy. 
The  results  show  a  relative  stability  of  the  ordered  alloy  for  which  the  lattice  constant  and  the 
bulk  modulus  are  predicted  to  be  6.20  A  and  53  GPa,  respectively.  Analysis  of  band  structure 
and  density  of  states  shows  the  alloy  to  be  direct-gap  semiconductor  with  the  value  of  4.44  for 
the  index  of  refraction  at  6889  nm. 
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First  Principles  Study  of  CdGeAso. 

1  Introduction 

Since  the  first  CdGeAso  crystals  were  grown®  in  the  early  70s,  their  bulk 
properties  have  attracted  much  attention.  The  photoluminescence  studies 
have  revealed  the  detailed  information  about  the  nature  of  its  band  gap.^ 
The  electrical  and  Hall  effect  meiisnrenients  have  determined  the  nature  of 
charge  carriers  in  this  material.®'®  The  presence  of  at  lejist  two  native  accep¬ 
tors  luus  b(!en  observed  in  experiments  utilizing  radiation  damage,*®  mag¬ 
netic  re.sonan<:e"  and  thermal  admittance  spectroscopy '*  techniciues.  Bulk 
tiiodulus  and  Debye  temperature  were  also  obtained  from  mcivsurements  of 
the  ultnusonic  wave  velocities.**  On  the  theonjtical  front,  empirical  pseu¬ 
dopotential  calculations  have  focused  only  on  its  band  structure. '■’•  *'*  Al¬ 
though  calculations  based  on  first-principles  methods  have  been  performed 
on  several  members  of  the  chalcopyrite  grotip,*^’*®  we  are  not  aware  of  any 
such  calculations  on  CdGeAs2. 

In  thisLsection,  we  report  the  results  of  calculations  based  on  the  density 
functional  theory  on  CdGeAs2  with  the  aim  (i)  to  provide  a  detailed  infor¬ 
mation  about  its  electronic  structure,  (ii)  to  determine  its  thermodynamic 
properties  within  a  Debye-like  model,  and  (iii)  to  combine  knowledge  of  elec¬ 
tronic  structure  and  thermodynamics  to  determine  pressure  dependence  of 
its  band  structure.  In  the  following  section,  computational  details  will  be 
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•riven.  Results  will  be  disccused  in  Sec.  IV. 

2  Computational  Details 

Electronic  structure  calculations  of  CdGeAso  are  based  on  the  density  func¬ 
tional  theory  embedded  in  the  CRYSTAL95  program^®  which  has  been  suc¬ 
cessful  in  describing  bulk  and  surface  properties  of  covalent  materials  such 
as  Si  and  GaN.'*^’"*  VVe  perform  calculations  in  the  framework  of  the  local 
density  approximation'-^'--  (LDA)  to  the  density-functional  theory  (DPT). 
In  some  cases,  gradient-corrected  functionals  (GG.A.),  Beclce'^'  for  exchange 
and  Perdew  and  VV.ing-’’' for  corndation,  an?  (iinployed  to  .see  the  effect  of 
tli(!.se  cornictions  on  the  calcidat(;d  structural  properties. 

A  lint;ar  combination  of  Gaussian  orbitals  are  used  to  generate  a  localized 
atomic  Ijiusis  from  which  Bloch  functions  are  constructed  by  a  further  linear 
combination  with  plane-wave  phase  factors.  To  reduce  the  computational 
co.st  in  the  present  work,  we  represent  the  core  electrons  by  Steven’s  effective 
core  potentials(ECPs)^"’^®  and  describe  the  respective  valence  electrons  by 
the  corresponding  Gaussian  basis  set.  The  core  electrons  for  Cd,  Ge  and 
As  are  Is^,  2s^,  2p*’,  3s^,  3p®  and  3d^°.  It  has  been  shown  that  the  ECPs 
(such  as  for  Zn  and  Ge  used  in  this  work)  derived  from  the  Hartree-Fock 
calculations  provide  satisfactory  description  of  the  ground  state  properties 
in  the  DPT  calculations.^®  The  Gaussian  basis  set  consists  of  four  s,  four  p 
and  three  d-type  shells  for  Cd  (i.e.  a  4121/4121/311  set),  two  s,  two  p  and 
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a  d-type  shell  for  Ge  (i.e.  a  41/41/1  set),  and  two  s,  two  p  and  a  d-type 
shell  for  As  (i.e.  a  41/41/1  set).  It  is  to  be  noted  here  that  the  outermost 
Ad  electrons  of  Cd  are  considered  explicitly  here  as  valence  electrons  and 
are  not  included  in  the  ECPs.  Furthermore,  the  outer  exponents  in  each 
of  the  birsis  set  are  reoptimized  to  yield  the  minimal  total  energy  at  the 
experimental  geometry  of  CdGeAso. 

The  initial  configuration  for  the  geometry  optimization  is  taken  from 
the  X-ray  diffraction  data.  At  ambient  conditions,  CdGeAs)  crystallizes 
in  the  chalcopyrite  phase^*^'^^  with  a  space  group  of  I-42D.  It  can  be  con¬ 
structed"  by  first  considering  a  superlattice  of  the  (cubic)  zincblende  phase 
with  »:/a=2,  replacing  «!ach  half  of  cations  by  Cd  ami  Ge  ions  respectively 
and  finally  introilncing  a  small  distortion  along  z-.ixis  leading  to  c/a  =  1.88!). 
Fnrtht’nnore,  the  As  atoms  are  placed  in  the  lattice  in  the  special  positions 
giv(;n  by  tin;  int(5rnal  parameter,  tt,  which  is  related  to  the  tetragonal  dis¬ 
tortion  (i.e.  2-(c/a))  in  the  lattice.  Geometry  optimization  calcidations  in 
the  LDA  approximation  therefore  involve  the  following  steps;  (i)  for  .several 
fixed  values  of  the  crystallographic  unit  cell  volume,  ranging  approximately 
from  0.8  to,  1.08  times  the  experimental  one,  the  lattice  constant  a  and  the 
internal  parameter  u  is  optimized  with  accuracies  better  than  10“^  A  and 
10“*,  respectively;  (ii)  the  calculated  potential  energy  surface  (i.e.  total 
energy  vs  volume)  is  then  used  to  obtain  the  optimized  value  of  the  unit 
cell  volume;  and  (iii)  the  optimization  of  a  and  u  is  again  performed  at  this 
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volume  to  get  the  values  of  equilibrium  structural  parameters.  In  the  GGA 
calculations,  geometry  optimization  is  achieved  by  a  less  -exhaustive  ap¬ 
proach  where  series  of  one  dimensional  searches  are  performed  varying  each 
structural  parameters  (i.e.  a,  c  and  u)  one  by  one  on  each  step,  repeating 
the  steps  until  convergence  of  0.01  A  in  the  parameter  value  is  achieved.  In 
all  of  these  calculations,  the  tolerance  on  the  total  energy  convergence  in  the 
iterative  solution  of  the  Kohn-Sham  equations  is  set  to  10“®  Hartree  and  a 
grill  of  50  Jt-points  was  used  in  the  irrudicible  Brillouin  Zone  for  integration 
in  the  reciprocal  space.^* 

3  Structural  and  Electronic  Properties 

The  calculated  eiiuilibrium  values  of  the  structural  parameters,  namely  a, 
c/a  and  u  for  CdGeA.s-j  are  5.014  A,  1.015  and  0.205  in  the  LDA  approxima¬ 
tion  and  5.04  A,  1.80  and  0.265  in  the  GGA  approximation  r(!spectivcly.  On 
the  other  hand,  the  experimental  values  of  the  lattice  parameters,^*  o  and  c, 
are  5.9432  A  and  11.2163  A  respectively  at  298  K  leading  to  the  c/a  ratio 
of  1.887.  Comparison  between  the  GGA  and  experimental  values  of  the  lat- 
tice  constants  shows  an  excellent  agreement.  The  calculated  LDA  unit  cell 
volume  is  also  well  within  1%  of  the  experimental  value.  However,  the  LDA 
calculations  yield  slight  underestimation  of  a  and  overestimation  the  cj a  ra¬ 
tio  as  expected.  Both  LDA  and  GGA  values  of  the  internal  parameter  u  are 
lower  than  the  experimental  value  of  0.2785.  This  is  in  accordance  with  the 
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results  of  the  LDA  calculations  on  ZnGePo,  CdSnP2  and  CdSnAs2  where 
the  calculated  u  is  consistently  lower  than  the  corresponding  experimental 
value.  It  is  to  be  noted  here  that  the  X-ray  diffraction  studies  obtain  the 
value  of  the  internal  parameter,  u,  only  in  an  indirect  way.^ 

The  calculated  interatomic  distance  between  Cd-As  and  Ge-As  is  2.579 
and  2.474  A  at  the  LDA  level  and  2.581  and  2.477  A  at  the  GGA  level, 
respectively.  It  has  been  suggested  that  the  tetrahedral  angle  in  the  chal- 
copyrite  lattice  is  generally  conserved  by  the  group-IV  atom  (i.e.  Ge).  This 
results  in  a  larger  separation  between  the  group  II  ami  the  group  V  atoms 
relative  to  the  group  IV  and  the  group.  V  atoms  in  the  lattice.  This  is  what 
hius  been  predicted  by  our  calculations  in  CdGeA-Sj.  A.ssuming  the  con- 
.servation  of  tetrahe<lral  bond  radii,  .Jalfe  and  Zungor*’'  hav<!  reported  the 
structural  coordinates  of  chalcopyritc^s.  In  this  schenu;,  the  calcidatcd  Cd- 
As  and  Ge-As  distance  comes  out  to  be  2.G3  and  2.45  A  respectively.  VVe 
note  here  that  the  LDA  calculations'*'*  in  CdSnAs2  yield  the  Cd-As  distance 
to  be  2.6 11  A. 

Fig.  1  displays  the  upper  valence  and  lower  conduction  band  structure 
of  CdGeAs2  calculated  at  the  LDA  equilibrium  geometry.  Total  and  partial 
density  of  states  of  the  valence  band  are  shown  in  Fig.  2.  Accordingly, 
the  upper  valence  band  is  mainly  formed  by  the  As  p  states  with  very  small 
admixture  of  the  Ge  p  states  and  the  Cd  p  states  (Fig.  2).  This  band 
has  a  width  of  about  5.5  eV.  In  the  absence  of  the  spin-orbit  interaction 
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terms  in  the  present  work,  the  top  of  the  valence  band  consists  of  two  levels 
having  symmetry  Tsu  and  r.i„;  the  former  being  double  degenerate.  Both 
levels  originate  from  the  same  triple  degenerate  F 15  level  of  the  sphalerite 
structure.  The  difference  between  the  Fsu  and  F-m  levels  is  referred  to  as. 
the  crystal-field  splitting  (Ac/)  which  is  calculated  to  be  0.33  eV.  The  cor¬ 
responding  experimental  value  is  reported  to  be  0.21  eV.' 

At  about  9  eV  below  the  top  of  the  valence  band  (Fig  2),  a  relatively 
narrow  Cd  Ad  band  overlaps  with  the  band  representing  the  bonding  be¬ 
tween  the  Ge  and  As  in  the  lattice.  Interestingly,  the  LDA  calculations 
on  CdSnPj  reported  the.  Cd  Ad  band  to  be  at  about  -8.5  eV  relative  to  the 
valinice  band  maximum.  Kowcwer,  we  note  here  that  the  DFT  calculations 
in  tin?  absence  of  .s(!lf-interaction  ccirrection  tend  to  underestimati!  the  loca¬ 
tion  of  the  d  band  with  reference  tb  the  top  of  the  valimcf?  IkIikI.  Fig.  2 
al.so  shows  a  band  at  about  -I I  eV  attributed  mainly  to  the  As  h  orbitals 
which  appcjar  to  be  hybridizctl  with  the  C«l  and  Ge  orbitals.  Relative  to 
the  valence  band  maximum,  we  have  collected  the  energies  of  the  electronic 
states  at  F,  T  and  tf  k-points^  that  correspond  to  maxima  and  minima  of 
the  valenCb  and  conduction  bands  in  Table  1. 

In  contrast  to  a  very  reasonable  description  of  valence  states,  one-electron 
solutions  of  Kohn-Sham  equations  corresponding  to  the  states  in  the  con¬ 
duction  band  are  known  to  be  much  less  reliable.  The  most  evident  and 
well-known  discrepancy  is  a  severe  underestimation  of  the  forbidden  band 
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gap.  For  example,  the  calculated  band  gap  in  the  present  work  for  CdGeASa 
is  0.16  eV  cis  compared  to  the  experimental  value  of  0.61  eV.  This  discrep¬ 
ancy  may  be  corrected  by  employing  a  semiempirical  correction^"*  which 
simply  shifts  states  in  the  conduction  band  using  the  following  equation  : 

9eV 


A  = 


(1) 


where  £oo  is  the  high-fretiuency  dielectric  constant  of  the  material.  For 
CdGeASo,  £oo  is  about  ILO*'--*"*  leading  to  the  value  of  the  gap  to  0.98  eV 
(Table  2). 

The  lowest  gap  between  top  of  the  valence  band  to  bottom  of  the  conduc¬ 
tion  band  is  found  to  lx;  at  T.  'rims,  the  present  DF'C  calctdations  predicts 
the  gap  to  be  direct  (r.|„-ric)  in  agreement  with  expc^riinental  and  empiri¬ 
cal  p.seudopotential  studies.*'*  The  ordering  of  states  in  the  conduction  band 
(i.(!.  Fic  <  F.i,;  <  r^c)  is  also  correctly  predicted.  We  note  here  that  exper¬ 
imental  assignments  arc  only  tentative  because  transitions  to  the  last  two 
states  arc  only  weakly  allowed. 


4  Thermodynamic  properties 

To  determine  thermodynamic  properties  of  CdGeAs2,  we  have  used  a  quasi¬ 
harmonic  Debye-like  model  in  which  the  Debye  temperature,  ©,  depends 
only  on  the  volume  of  the  crystal.  It  is  a  non-empirical  model  that  only 
needs  the  knowledge  of  the  potential  energy  surface  (i.e.  a  set  of  computed 
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{V,  Eiaiti^))  points)  to  generate  the  equation  of  state  (EOS)  and  a  number 
of  related  quantities  of  a  given  material.  It  is  relevant  to  notice  that  whereas 
the  key  parameters  from  the  experimental  point  of  view  are  pressure  (p)  and 
temperature  (T),  the  appropriate  variable  to  control  the  computation  is 
volume  (K).  The  computational  strategy  therefore  includes  optimization  of 
the  lattice  constant  a  and  the  internal  parameter  u  for  several  fixed  values  of 
the  unit  cell  volume  to  obtain  the  energy  surface,  (see  Sec.  2).  Fig.  3  shows 
such  energy  surface  for  CdGeAs.)  calculated  using  the  LDA  approximation. 

Assuming  the  isotropic  conditions,  ©  can  be  given  by^^ 

where  /i  is  the  reduee<l  Planck  constant,  ku  is  the  Boltzmann  constant,  M 
the  molecular  tn;uss  of  tlie  compound,  r  the  number  of  atoms  per  molecular 
unit,  Z7.V  the  adiabatic  I)ulk  modulus  of  the  crystal,  and  cr  the  Pois.son  ratio 
which  is  taken  to  be  0.25,  the  value  of  the  Cauchy  solid.'^''’  The  value  of 
/((T  =0.25)  then  becomes  0.85995.  Note  that  the  explicit  expression  for 
/(or)  in  Eq.  2  is’’'' 


In  Eq.  2,  Bs  depends  on  V  and  T.  In  order  to  balance  computational 
demand  and  accuracy,  we  further  assume 
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where  ZJstaiic  >s  the  static  bulk  modulus. 

0  is  now  simply  a  function  of  V.  It  is  also  evident  from  Eqs.  2  and  4  that 
there  is  no  explicit  consideration  of  the  structural  parameters  in  the  volume 
derivatives,  since  we  restrict  ourselves  to  hydrostatic  conditions.  This  rele¬ 
vant  feature  of  the  model  makes  it  fully  independent  of  any  particular  crystal 
structure.  Even  if  the  crystal  energy  depends  on  many  internal  parameters 
and  cell  constants,  a  set  of  pairs  (Eiatt,  V)  is  sufficient  to  employ  the  model. 
A  more  detailed  description  of  the  model  is  given  in  Ref.  (34,35). 

Our  computed  values  of  isothermal  and  adiabatic  bulk  moduli  (Brip  = 
0)  =  /?()  aiul  Ds)  itiul  their  pressure  and  temperature  derivatives  are  given 
in  'lable  3  along  with  the  availablt;  (!xp(;rimental  data.  Fig.  4  shows  the 
g(merated  etpiation  of  state  tising  the  normalized  V(p)/Vi)  values,  where  l^o 
is  the  corresponding  etpdlibrium  volume  at  zero  pressure.  It  shotdd  be  noted 
here  that  Hailing  at  «/*'*  used  the  experimental  bulk  modulus  and  its  first 

s 

derivative  obtained  .at  pressures  up  to  0.1  GPa  from  ultr.xsotdc  me.osurements 
.as  p.ar.ametcrs  and  did  not  fit  the  equation  of  state  to  the  cxpcriment.al  curve. 
The  higher  experimental  value  of  69.7  GPa  of  Bq  rel.ative  to  the  calculated 
one  of  57.6  GPa  makes  the  simulated  crystal  to  be  more  compressible,  as 
illustrated  in  Fig.  4. 

In  the  quasi-harmonic  model  used  here,  the  Debye  temperature,  0Oi 
comes  out  to  be  273K  at  the  static  equilibrium  and  269K  at  room  temper¬ 
ature.  This  vzuriation  is  due  to  change  of  equilibrium  volume  with  temper- 
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ature.  On  the  other  hand,  the  specific  heat  mecisurernents^^  between  2K 
and  30K  determine  Qq  to  be  230K  whereas  the  ultrasonic  wave  velocity 
experiments^  yields  the  value  of  257K  for  ©o- 

We  note  here  that  thermodynamic  properties  that  are  specially  sensi¬ 
tive  to  computational  details  are  those  related  to  the  volume  dependence 
of  0,  since  they  involve  third  derivatives  of  the  static  lattice  energy.  The 
most  important  factor  in  computing  this  set  of  properties  is  the  Griineisen 
constant,  7.  Our  calculations  yield  7  to  be  2.2(12  at  300K  lor  CdGeAs). 
For  the  volume  thermal  expansion  coefficient,  ct,  the  calculatetl  value  is  30.9 
X  10"*’  K~s  at  3()0K,  and  33.8  x  l()~®  K~‘  at  dOOK  :us  compared  to  the 
e.xpirrimenlal  value  ot  18  x  lO'*'  K“’. 

The  optical  transitions  in  chalcopyrites  are  .sensitive  to  applied  hy<lro- 
static  pressure  Icsuling  to  dependence  of  the  pre.s.snr(r coefficients  on  the 
nature  of  the  band  gap.  For  example,  it  h.as  been  found'”s  that  the  p.sendo- 
diroct  gap  chalcopyrite.s,‘^ss  such  as  ZnGeP2,  have  a  small  positive  or  negative 
slope  of  the  absorption  edge,  while  CdGeAs^  and  other  direct  gap  chalcopy¬ 
rites  have  a  large  positive  slope.  Although  the  LDA  calctdations  do  not 
provide  quantitative  predictions  of  the  band  gap  very  well,  they  do  provide 
reliable  predictions  of  the  change  in  gap  under  the  hydrostatic  conditions. 
For  CdGeAs2,  we  find  that  both  valence  and  conduction  levels  at  T,  iV  and 
r  k-points  in  the  Brillouin  zone  show  a  linear  vziriation  with  applied  pressure 
up  to  10  GPa.  The  calculated  pressure  coefficients'^®  of  these  levels  are  given 
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in  the  Table  4  which  yield  the  coefficient  of  the  minimum-energy  direct  gap 
(r.iu-Fic)  to  be  9.5  x  10"®  eV/bar.  It  is  in  very  good  agreement  with  the 
experimental  value'^  of  9.3  x  10"®  eV/bar.  It  is  expected  that  knowledge 
of  the  pressure  coefficients  for  different  levels  would  assist  experimentalists 
in  identifying  various  optical  transitions  in  CdGeAso. 

In  summary,  we  have  calculated  lattice  constants,  etpiation  of  state, 
Debye  temperature,  Griineisen  constant  aiul  electronic  band  structure  of 
CdGeAso  in  chalcopyrite  structure  which  are  in  very  good  agi’eement  with 
the  availble  experimental  data.  The  calculatetl  results  have  also  provided 
details  of  the  valence  and  conduction  bands  predicting  their  variation  with 
pre.ssure. 
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Table  1:  Valence  band  energies  in  oV  at  F,  T,  N  k-points  in  the  Brillouin 


zone  of  GdGeAs-j. 


Level 

As  p  band 

maxima 

r4u 

0.00 

Tsu 

-0.33 

7;$u  +  T.lu 

-1.26 

Niu 

-0.72 

As  p  baud 

minima 

r.,„ 

,4.49 

T.ly  +  Tju 

-3.86 

/v,„  . 

-.'j.OO 

(Ji!-As  band 

V.U. 

-7.0.1 

r*2r/ 

-9.44 

Tr,„ 

-7..1() 

/v,„ 

-7..'32 

Cd  fl  band 

r  rnajc 

-8.24 

r  min 

-8.8.5 

Trnax 

-8.40 

Tmin 

-8.72 

^max 

-8.50 

^min 

-8.90 

As  s  band 

r.5„ 

-11.04 

Tsv 

-11.09 

ri„ 

-13.10 

Tlu  "b  Tsu 
T^y 
Nly 

Nly 


6^ 


-11.09 

-12.34 

-11.23 

-12.19 
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Table  2:  The  minimum-energy  direct  and  pseudo-direct  band  gaps  and  crys¬ 
tal  field  splittings  (in  eV)  for  CdGeAso 


DFT  gap 

Corrected  gap 

Experiment 

^direct 

0.16 

0.98 

0.61“ 

^pdtuidu^dircct 

1.16 

1.98 

■ — 

Acf 

-0.33 

— 

-0.21 

a  Ref.  (7] 


Table  3:  Bulk  moduli  and  related  properties  of  CdGeAso 


BoCGPa) 

5s{GPa) 

B'o 

f??(GPa-') 

This  work  :  OK 

57.03 

57.03 

4.79 

-0.32 

Tlus  work  :  300K 

53.21 

54.43 

5.47 

-0.46 

Experiment":  300K 

69.7 

: - 

6.18 

. - 

Table  4:  Pressure  coefficients  in  (10  **  eV/bar)  of  valence  and  conduction 
levels  of  CdGeAS). 


Level 

Tr,u 

o.y 

r,„ 

G.5 

r;,.: 

-2.0 

Vic 

1() 

T'Iu+'Im 

-0.87 

2.2 

2.6 

iV.c 

8.4 

Figure  Captions 

Fig.  I  :  Electronic  bcind  structure  of  CclGeAsj. 

Fig.  2  :  Total  atul  Projected  density  of  states  of  CdGeAs  >. 

Fig.  3  :  Total  Energy  vs  volume  of  CdGeAso  in  chalcopyrite  phase 
Fig.  4  :  Etiuation  of  State  of  CdGeAs>. 
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Atomistic  Modeling  of  Native  Defects  in  CdGeAs2 


I.  INTRODUCTION  : 

Cadmium  germanium  arsenide  (CdGeAs2)  is  an  important  material  due  to  its 
suitability  for  nonlinear  optical  applications  in  the  infrared  region  since  it  has  the  highest 
nonlinear  optical  coefficient,  236  pm/V,  known  for  a  phase  matchable  compound 
semiconductor.  Additionally,  it  has  a  wide  transparency  region  extending  across  the 
infrared  from  2.4  to  18  The  availability  of  large  single  crystals^*  ofCdGeAs2  has 

led  to  renewed  efforts  to  utilize  this  material  for  frequency  doubling  (i.e.,  for  second 
harmonic  generation)  carbon  dioxide  laser  emission  lines  from  9.3  to  10.6  (xm.  The 
transparency  of  those  crystals'*,  although  superior  to  that  of  earlier  crystals'  is  still  limited 
by  the  photo-ionization  of  native  acceptors. 

Nominally  undopcd  Bridgman  grown  single  crystals  of  CdGeA.S2  are  found  to  be 
p-type.  Hall  effect  measurements^-^’  on  the  best  available  crystals  indicate  the  pre.sence  of  a 
single  acceptor  with  an  activation  energy  of  0. 15  eV.  On  the  other  hand,  Brudni  et.  al."^, 
Halliburton  et.  al.^  and  Smith  et.  al.^  utilizing  re.spectively  radiation  damage  techniques, 
EPR  and  thermal  admittance  spectro.scopy  measurements  observe  that  there  can  be  at  least 
two  native  acceptors  present  in  as-grown  CdGeAs2 .  Smith  et.  al.^  have  found  that  there 
is  a  level  with  an  activation  energy  of  0.1 1-0.13  eV  in  agreement  with  Hall  measurement 
on  the  same  samples,  and  a  second  deeper  acceptor  with  an  activation  energy  of  0.346  eV. 
It  is  the  goal  of  this  work  to  determine  theoretically  which  native  defects  might  give  rise  to 
these  two  acceptor  centers  in  CdGeAs2. 

Theory  has  so  far  focused  on  calculating  intrinsic  properties  such  as  the  band 
structure'®-"  and  has  not  considered  the  chemistry  of  defects  in  this  material.  In  this 


paper.  we  perform  such  a  task  considering  vacancies,  interstitials  and  antisites  in 
CdGeAs2.  We  use  atomistic  simulation  methods  based  on  the  shell  model  to  describe  the 
crystalline  lattice  of  CdGeAs2  and  will  calculate  formation  energies  of  native  defects  and 
binding  energies  of  the  acceptor  centers.  We  note  here  that  this  approach  has  been  very 
successful  in  yielding  reliable  defect  energies  in  a  wide  variety  of  materials  including 
sixfold-  and  fourfold-coordinated  structures. ^2  Recent  applications  of  the  shell  model 
include  structural  relaxation  around  dopants  in  sapphire '3,  high-pressure  phase  transition  in 
GaN  and  defect  energetics  in  ZnGeP2 . 

11.  PERFECT  LATTICE 

In  the  present  atomistic  description  of  a  crystal,  the  lattice  is  considered  to  consist  of 
ions  interacting  via  a  combination  of  electrostatics  and  interatomic  potentials,  while 
polarization  is  included  by  means  of  the  shell  model.  The  total  energy  of  a  cry.stal  is  then 
taken  to  be  a  sum  of  two-  and  three-body  interactions  in  the  lattice.  The  (two-body) 
pairwise  interaction  term  consists  of  the  long-range  Coulombic  part  and  the  short-range 
repulsive  part  and  is  given  by 

Eij  =  (QiQj/Rij)-H  VsR(Rij).  (I) 

where  Q's  are  charges  and  Rjj  is  the  separation  between  ions.  Vsr  (Rjj )  repre.sents  the 
short-range  interaction  between  ions  and  is  given  by  an  analytical  expre.ssion  of  the  Born- 
Mayer  form :  ^ 

^SR  (^ij )  ~  Aij  exp  ( -Rij  /  pij),  (2) 

where  A  and  p  are  the  parameters  obtained  generally  by  empirical  fitting  methods. 

Three-body  interaction  terms  in  our  model  are  represented  by  the  Axilrod-Teller*^ 
potentials: 

Ejjk  =  kijk  (1+3  cos  0i  cos  0j  cos  0ic )  /  Ri^  Rj^  Rk^.  (3) 
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where  kjjk  is  a  coefficient,  dj  is  the  ith  angle  and  Ri  is  the  side  of  the  triangle  formed  by 
ions/. yam/ in  the  lattice. 

To  describe  the  polarizabilty,  we  employ  the  shell  model  description  of  the  lattice 
atoms.  In  the  shell-model,  each  atom  consists  of  a  core  of  charge  X,  and  a  shell  of  charge 
Y,  such  that  the  total  charge  is  the  sum  of  the  core  and  shell  charges.  The  polarization  is 
then  simulated  by  the  displacement  of  a  shell  from  a  core,  the  two  being  connected  by  a 
harmonic  spring  with  a  force  constant  K. 

CdGeAs2  crystallize  in  the  chalcopyrite  phase' with  a  tetragonal  symmetry 
group  of  I42d.  Chalcopyrites  can  be  constructed-'^  by  first  considering  a  superlattice  of  the 
(cubic)  zincblende  phase  with  c/a=2,  replacing  each  half  of  cations  by  Cd  and  Ge  ions 
respectively  and  finally  introducing  a  small  distortion  along  z-axis  leading  to  c/a  =  1.889. 
As  shown  in  Fig.  1.  each  atom  is  tetrahedrally  coordinated  in  the  lattice;  cadmium  or 
uermaniiim  cations  have  4  near-neighbor  arsenic  anions  while  arsenic  anions  have  2 
cadmium  and  2  germanium  cations  as  near-neighbors.  In  the  lattice,  cations  occupy  special 
positions  labeled  as  a  and  h  with  no  positional  degrees  of  freedom,  while  anions  arc 
placed  in  the  cl  positions  given  by  the  internal  parameter,  x  .  The  value  of.r  varies  from 
crystal  to  crystal  and  is  related  to  tetragonal  distortion  (i.c.  2  -  c/a)  in  the  lattice.  For 
CdGeAsi  the  lattice  parameters '  9,  a  and  c  are  5.9432  A  and  1 1.2163  A  respectively  at 
298  K  and  the  parameter  .r  is  0.2785  leading  to  the  interatomic  distances  of  2.629  A  for 
Cd-As  and  2.430.  A  for  Ge-As. 

In  the  chalcopyrite  lattice,  the  tetrahedral  coordination  of  atoms  suggests  that  the 
covalent  bonding  (with  sp^-hybrid  bonds)  predominates.  On  the  other  hand,  the 
composition  of  the  cation  sublattice  indicates  a  presence  of  the  ionic  character  in  the 
chemical  bonding.  In  CdGeAs2  a  relatively  large  magnitude  of  tetragonal  compression 
(=6%)  confirms  that  the  chemical  bonding  is  mixed.  Based  on  the  electronegativities  of  the 
constituent  atoms,  the  Ge-As  bond  is  expected  to  be  less  ionic  than  the  Cd-As  bond.  Our 
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atomistic  model  therefore  do  not  assume  either  fully  covalent  (i.e.  Cd-‘Ge‘^As2'‘^)  or  fully 
ionic  (i.e.  Cd-‘'’Ge'^’*'As2^')  description  of  the  lattice,  but  uses  the  empirical  fitting  method 
to  determine  the  charges  associated  with  Cd,  Ge  and  As.  Hence,  the  mixed  (ionic-covalent) 
chemical  bonding  is  taken  into  account  by  the  use  of  the  derived  atomic  charges  along  with 
the  three-body  terms  included  in  the  atomistic  description  of  the  material. 

In  our  model,  the  two-body  interaction  term  (Eq.  2)  therefore  requires  the 
determination  of  Qcd.  Qce.  Qas  iind  Vsr  (i.  e.  parameters  A  and  p )  for  the  Cd-As,  Ge-As 
and  As-As  interactions  in  CdGeAs2.  Note  that  the  short-range  interactions  between  Cd-Cd 
and  Ge-Ge  are  ignored  as  they  become  very  small  for  large  separations  (>  4.0  A)  in  the 
lattice.  We  also  need  to  determine  the  coefficient  k^sGlAx  It)*"  the  three-body 

terms  (Eq.  3)  along  with  the  shell-model  parameters  Yas  and  kAs-  Both  cations  are 
considered  as  rigid  ions  in  the  lattice.  All  of  these  potential  parameters  are  obtained  by  the 
empirical  fitting  method-’  .which  u.ses  the  experimentally  known  crystal  properties  .such  as 
crystal  structure clastic-- and  dielectric constants  of  CdGoAs2.  Fitting  and  all 
calculations  were  performed  using  the  program  GULP.-'^ 

Table  I  lists  the  model  parameters  representing  the  interatomic  interactions  in  the 
lattice.  The  calculated  crystal  properties  arc  compared  with  the  experimental  data  in  Table 

II.  Accordingly,  the  potential  model  reproduces  the  lattice  structure  very  well.  The  overall 
good  agreement  between  the  calculated  and  experimental  properties  for  the  perfect  lattice 
indicates  the  reliability  of  the  derived  interatomic  potential  set  for  CdGeAs2. 

III.  NATIVE  DEFECTS 

Defect  energies  of  several  plausible  types  of  ionic  and  electronic  native  acceptor 
defects  have  been  calculated  using  the  Mott-Littleton  methodology ^5.26  in  which  the  lattice 
is  divided  into  a  series  of  different  regions  around  the  defect  by  concentric  spherical 
boundaries.  Immediately  surrounding  the  defect  is  region  1  in  which  all  ions  are  treated 
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explicitly  and  fully  optimised.  Beyond  thi.s  is  region  2a  in  which  all  atoms  arc  still  explicitly 
considered,  but  the  relaxation  effects  are  much  smaller  and  can  be  treated  more 
approximately.  The  validity  of  the  successive  approximations  made  in  the  above 
multiregion  method  improves  as  the  radii  of  the  regions  increase.  Therefore  it  is  important 
to  check  the  convergence  of  the  results  with  region  size.  In  the  present  work  a  region  1 
containing  approximately  350  atoms  was  found  to  be  sufficient  to  converge  the  absolute 
defect  energy  to  better  than  0.0 1  eV,  though  relative  energies  will  be  far  more  converged 
than  this.  We  also  note  here  that  provide  the  defects  are  charge  neutral  overall  and  the 


quality  of  fitting  to  the  structural  properties  are  the  same,  then  to  the  first  order  the  defect 
energies  should  not  depend  too  strongly  on  the  charges  used  in  the  model.  If  defect- 
induced  displacements  in  the  lattice  are  large  then  the  particular  shape  of  the  anharmonic 
region  of  the  potential  energy  curve  will  become  more  important  thereby  yielding  a  strong 
dependence  of  defect  energies  on  the  charges.  However,  given  that  we  know  the  ions  in 
CdGcAsT  do  not  have  integral  charges,  the  u.se  of  partial  charges  is  therefore  expected  to 
give  a  better  account  of  this  anharmonic  region  in  defect  calculations. 

The  Schottky  defects  in  the  lattice  are  formed  by  moving  the  constituent  ions  to  the 


surface  from  their  bulk  sites  whereas  the  Frenkel  defects  are  pairs  of  the  vacancy  and 
interstitial  of  the  same  type  of  ion.  In  calculations  an  interstitial  was  simulated  as  an 
addition  of  atom  to  the  empty  interstitial  position  in  the  lattice.  The  Schottky,  Frenkel  and 
antisite  formation  energies  (Table  III)  were  obtained  from  defect  energy  calculations  of 
vacancies,  interstitials  and  antisites  in  the  lattice.  In  CdGeAs2  the  Schottky  defect  is 
(Vcd+VGe+2VAs).  the  Frenkel  defect  pairs  are  (Vcd+Cdi),  (Vce+Gei)  and  (Vas+Asi)  and 
the  antisite  pairs  are  (Cdoe+Gecd)*  (GdAs+Ascd)  ttnd  (GcAs+ASGe)-  As  shown  in  Table 
ITT,  the  lowest  formation  energy  comes  out  to  be  for  the  (CdGe+Cecd)  antisite  pair  which  is 
followed  by  the  Schottky  and  Frenkel  pairs  of  Cd  and  Ge.  The  large  formation  energies  for 
the  antisite  pairs,  (CdAs+Ascd)  and  (Gcas+Asgc).  would  seem  to  preclude  their  occurrence 
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as  intrinsic  point  (Jcfccts  in  CdGcAsi.  This  ordering  of  formation  energy  is  what  we 
expected  since  the  magnitude  of  formation  energy  depends  mainly  on  the  extent  of 
distortion  introduced  by  individual  defects  in  the  lattice.  In  the  present  case,  antisite 
disorder  in  the  cation  sublattice  introduces  least  distortion  in  the  lattice  compared  to  that 
introduced  by  either  vacancies,  interstitials  and  antisites  involving  both  cation  and  anion 
sublattices. We  note  here  that  (tetrahedral)  covalent  radii  of  Cd,  Ge  and  As  ions  are 

1.405,  1.225  and  1.225  A  respectively  27. 

The  magnitude  of  the  formation  energy  for  the  (Cdce+Gecd)  antisites  is  small 
(‘=‘0.3  eV)  suggesting  that  appreciable  disorder  would  occur  in  the  cation  sublattice  at 
higher  temperatures.  This  is  in  accord  with  what  has  been  reported  for  CdGeAs2  which  has 
a  high-temperature  disordered  zincblende  phase.  In  this  phase,  Cd  and  Ge  occupy  cation 
sites  randomly  in  contrast  to  the  ordered  chalcopyrite  phase  where  Cd  and  Ge  occupy 
alternating  cation  sites.  Upon  cooling  to  room  temperature,  the  lattice  i.s"  then  expected  to 
retain  some  of  the  disorder  in  the  cation  sublattice. 

IV.  NATIVE  ACCEPTORS  : 

A  native  acceptor  center  in  CdGeA.S2  may  involve  either  cation  vacancy  (Vcj  or 
Vgc)  or  anlisitc  (Cdcc)  ba.sed  on  the  effective  charge  considerations  and  can  stabilize  a  hole 
on  near-neighbor  A.y  ions  with  the  binding  energy  given  by  : 

AE  =  Ea  +  Eh  -  Eh+A  W 

where  Eh+a  refers  to  energy  of  the  hole-acceptor  complex.  Eh  is  energy  of  the  hole  and  Ea 
is  energy  of  the  acceptor  in  the  lattice.  For  the  bound  hole,  AE  is  taken  to  be  positive. 

The  hole-acceptor  complex  in  CdGeAs2  can  be  classified  according  to  the  degree  of 
delocalization  of  the  hole  states  in  the  lattice.  For  example,  a  complete  delocalization  refers 


a  a  detailed  list  of  individual  defect  energies  and  lattice  coordinates  can  be  obtained  from  the  authors 
(pandey@mtu.edu) 
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to  sharing  of  a  hole  by  all  of  the  four  A.v  neighbors  of  the  acceptor  (i.e.  4-center  case). 
Alternatively,  a  hole  ne.xt  to  the  acceptor  center  may  be  shared  by  a  pair  of  Ai’  neighbors 
(i.e.  2-center  case)  or  localized  on  one  of  the  Ai-  ions  (i.e.  I-center  case)  in  the  lattice.  The 
presence  of  hole-acceptor  comple,\e.s  in  as-grown  CdGeAs2  has  been  revealed  by  Hall- 
effect  measurements^-^’-^,  electron  paramagnetic  resonance  (EPR)^  ,  and  thermal 
admittance  spectroscopy.^  We  note  here  that  a  direct  comparison  can  be  made  between  AE 
and  the  acceptor  binding  energy  obtained  from  the  Hall-effect  measurements. 

Our  approach  to  calculation  of  AE  is  similar  to  that  used  in  an  earlier  study  of 

ZnOePi  where  we  assumed  that  the  presence  of  a  hole  in  the  lattice  reduces  the  (shell) 
charge  thereby  only  modifying  long-range  Coulombic  interactions  between  ions.  The  short- 
range  interaction  parameters  (Eq.  2)  are  taken  to  be  the  same  as  obtained  for  the  perfect 
lattice.  This  assumption  is  likely  to  overestimate  the  1-center  hole  binding  energy  slightly 
since  it  provides  li  more  accurate  treatment  of  the  delocalized  hole  (4-center  case)  relative  to 
the  localized  hole  ( l-center  ca.se).  Furthermore,  we  e.xpect.  to  a  first  appro.ximation, 
cancellation  of  the  hole  .self-interaction  energies  in  Eq.  4  for  the  4-center  ca.se  since  dE  is 
the  difference  of  two  terms  involving  the  delocalized  hole. 

Calculations  ba.sed  on  the  Mott-Littleton  methodology  (de.scribed  in  Sec.  Ill)  are 
performed  to  obtain  defect  energies  for  eq.  4  to  calculate  AE.  The.se  include,  for  example, 
antisite  (Ea),  delocalized  hole  (Eh)  and  antisite  next  to  the  delocalized  hole  (Eh+a)  the 
lattice.  All  defect  coordinates  are  fully  optimized  in  these  calculations.  The  calculated 
binding  energies  (AE)  are  given  in  Table  IV.  Accordingly,  a  striking  distinction  is  found  in 
the  binding  energies  of  acceptors  associated  with  either  vacancies  or  antisites  for  different 
hole  localization  regions.  The  hole  tends  to  be  localized  near  a  cation  vacancy,  i.e.,  the 
more  the  hole  is  localized,  the  larger  is  the  binding  energy.  But  the  trend  is  exactly  opposite 
for  antisites  where  a  larger  binding  energy  is  found  to  be  associated  with  the  delocalized 
hole.  It  is  well  known  that  the  localization  of  a  hole  in  the  lattice  is  a  result  of  interplay 
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between  lattice  distortion  and  polarization.  In  CdGeAsi.  the  distinct  behavior  of  hole 
localization  can  therefore  be  understood  in  terms  of  the  distortion  and  polarization 
introduced  by  these  defects  in  the  .surrounding  lattice.  Fore.xample,  a  Vcd  has  an  effective 
charge  of -1.8e  in  comparison  to  the  effective  charge  of  -0.5e  associated  with  Cdce-  The 
difference  in  effective  charges  along  with  the  fact  that  a  Cd  atom  is  much  larger  than  the  Ge 
atom  results  in  an  inward  relaxation  (=6%)  of  As  neighbors  for  Vcd  in  contrast  to  an 
outward  relaxation  (=5%)  of  As  neighbors  for  Cdce-  Based  on  simple  electrostatic 
arguments,  the  localized  hole  (i.e.  1-center  case)  is  therefore  expected  to  have  a  larger 
binding  energy  for  Vcd  than  that  for  Cdoe- 

As  shown  in  Table  IV  the  calculated  binding  energies  for  various  acceptor  centers 
vary  from  0.04  to  0.80  eV.  Note  that  the  energy  gap  for  this  material  is  only  0.57  eV  at 
room  temperature.  For  Hall  effect  measurements  on  the  best  available  crystals^*^,  a  value  of 
0. 15  eV  is  obtained  for  the  acceptor  binding  energy.  The  acceptor  center  revealed  in  the 


Hall-effect  measurements  is  therefore  likely  to  be  a.s.sociated  with  the  delocalized  hole  (i.e. 


4-center  ca.se)  in  the  vicinity  of  Cdcc  which  calculations  yield  the  binding  energy  of 


about  0. 13  eV.  The  acceptor  center  involving  Vcd  was  not  cho.sen  as  one  would  expect  the 
deepest  level.  0.49  eV,  to  be  its  populated  level  (Table  IV).  We  note  here  that  a  detailed 
electronic  structure  .study  of  the.se  defect-complexes  is  planned  to  confirm  the  proposed 
assignement  of  the  acceptor  centers  in  CdGeAs2. 

In  fact,  our  prediction  regarding  the  nature  of  the  dominant  acceptors  is  confirmed 
by  the  EPR  measurements.  Halliburton  et  al.^  have  proposed  that  the  dominant  EPR  signal 
may  be  due  to  the  acceptor  center  consisting  of  either  a  cation  vacancy  or  antisite  cation 
with  the  unpaired  spin  shared  by  the  four  neighboring  As  ions.  This  is  consistent  with  our 
selection  of  the  shallow  acceptor  as  the  cation  antisite,  Cdce-  Our  further  analysis  of  the 
calculated  results  find  that  the  near-neighbors  of  Vcd  relax  inward  by  0.08  A  in  contrast  to 
the  outward  relaxation  of  near-neighbors  of  Cdoe  of  about  0.  19  A  from  their  regular 


lattice  sites  (Table  V).  It  is  to  be  noted  here  that  the  magnitude  of  the  lattice  distortion 
introduced  by  acceptors  can  be  used  in  the  analysis  of  (future)  magnetic  resonance  studies 
to  ascertain  more  precisely  the  nature  of  the  EPR-active  center: 

Recent  capacitance-voltage  measurements  using  thermal  admittance  spectroscopy^ 
have  reported  the  presence  of  two  acceptor  levels  in  p-CdGeAsi  with  binding  energies  of 
0  1 1-0.13  eV  and  of  0.346  eV.  The  first  acceptor  level  is  the  same  one  revealed  by  Hall 
effect  measurements  yielding  binding  energies  of  0.10  -  0.12  eV  for  these  exact  same 
samples^  which  are  fairly  heavily  doped.  This  value  corresponds  to  the  delocalized  hole 
(i.e.  4-center  case)  bound  to  the  Cdce  tts  previously  discussed.  Based  on  comparison  of  the 
binding  energy  of  the  second  acceptor  with  the  calculated  binding  energies  (Table  V),  we 
predict  that  the  deeper  acceptor  is  associated  with  the  localized  hole  bound  to  Vc<i,  as  this 
center’s  deepest  level  will  be  populated.  The  other  possibilities  can  be  ruled  out  as  they 
have  binding  energies  which  e.xceed  the  band  gap  of  (TdGeAsi. 


5.  SUMMARY 

A  set  of  interatomic  potentials  consisting  of  two-  and  three-body  potential 
parameters  coupled  with  the  shell  model  description  of  the  lattice  is  developed  for  CdGeA.S2 
reproducing  its  structural,  elastic  and  dielectric  properties  .succe.ss fully.  Native  defect 
formation  energies  show  that  antisites  in  the  cation  sublattice  will  dominate  the  intrinsic 
disorder  in  this  material.  The  dominate  acceptor  center  controlling  the  optical  properties  of 
this  material  and  revealed  in  the  EPR  and  Hall-effect  measurements  is  likely  to  be  a  hole 
shared  by  four  neighboring  Aj  ions  in  the  vicinity  of  Cdce-  The  deeper  level  is  probably 
due  to  the  a  hole  localized  on  one  As  ion  and  bound  to  Vcd  • 

Although  ZnGeP2  and  CdGeAs2  belong  to  the  same  chalcopyrite  family,  the  nature 
of  dominant  acceptors  in  these  materials  is  predicted  to  be  different.  Defect  calculations 
corroborating  the  Hall  effect  and  magnetic  resonance  studies  find  that  the  zinc  vacancy,  not 
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the  zinc  antisite  (Zncc)*  >-'>  associated  with  the  dominant  acceptor  center  in  2nGeP2.  This  is 
not  the  case  in  CdGeAs2  where  a  cadmium  antisite  (Cdoc)  is  predicted  to  be  associated 
with  the  dominant  acceptor  in  the  lattice.  We  believe  that  this  difference  may  well  be  due  to 
defect-induced  lattice  distonion  which  plays  a  key  role  in  stabilizing  the  hole  states  in  the 
lattice.  Based  on  the  size  argument,  native  acceptors  associated  with  zinc  antisites  (Znoe) 
are  not  expected  to  introduce  significant  lattice  distortion  In  ZnGeP2  since  both  Zn  and  Ge 
have  the  same  tetrahedral  radius^^  of  1.23  A.  On  the  other  hand,  the  tetrahedral  radius  of 
Cd  is  about  1.41  A.  Therefore,  the  size  difference  of  about  15%  between  Cd  and  Ge  would 
be  expected  to  cause  significant  distortion  by  a  cadmium  antisite  (Cdco)  in  the  CdGeAs2 
lattice. 
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TABLE  I.  Interatomic  potential  parameters  of  Eqs.  2  and  3.  The  atomic  charges  are 
+  I.83e,  +2.29e  and  -2.06e  for  Cd,  Ge  and  As,  respectively.  For  As.  the  shell  charge  (Y) 
is -2. 22e  and  the  spring  constant  (K)  is  1.40  eVA'-. 


Aij  (eV) 

Pij  (A) 

kijk(eVA-9) 

Cd- As 

329.5 

0.4826 

Ge  -  As 

531.6 

0.4119 

As  -  As 

Cd-As-As 

32228.9 

0.3131 

272.4 

Ge-As-As 


-277.5 


TABLE  II.  Calculated  and  Experimental  bulk  properties  of  CdGeAsi . 


11 

11 

11 

11 

11 

11 

11 

11 

11 

11 

11 

11 

11 

11 

11 

11 

11 

It 

11 

This  work 

Experiment 

0 

Structure  a,  A 

5.951 

5.9432“ 

a 

c,  A 

11.221 

11.2163 

c/a 

1.886 

1.887 

XAs 

0.279 

0.2785 

Lattice  energy/unit  cell,  eV 

El  -134.6 

Elastic  constants  ( 10*  •  dyn  cm-2) 

Cm 

9.70 

9.45'^ 

Ci2 

4.47 

5.96  • 

C|3 

5.45 

5.97 

Cj3 

8.95 

8.34 

C44 

4.26 

4.21 

C66 

3.42 

4.08 

Dielectric  constants 

EO" 

14.75 

14.8*^ 

eoj3 

15.43 

15.4 

£00“ 

10.17 

10.06<‘ 

10.93 

11.04 

(a)  reference  19  (b)  reference  22 

(c)  reference  23 

(d)  reference  1 
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TABLE  III.  Native  defect  formation  energies  in  CdGeAsi . 


Formation  energy  (per  defect),  eV 


Schottky  pair  : 

(Vcd+Vce+2VAs) 

1.88 

Frenkel  pair  : 

(Vcd+Cdi) 

2.39 

(Vce+Gej) 

2.82 

(Vas+As-,) 

5.07 

Antisite  pair  : 

(Cdcc+Gecd) 

0.28 

(CdAs+A.sar) 

12.1 

(Gcas+Asgc) 

13.1 
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TABLE  IV.  The  calculated  binding  energies  (AE  of  eq.  4)  of  native  acceptors  in  CdGeAs2. 


Acceptor  center" 

1 -center 

2-center 

4-center 

Hole  near  Vcd 

0.49  eV 

0.26  eV 

0.15  eV 

Hole  near  Vce 

0.80 

0.67 

0.59 

Hole  near  Cdce 

0.04 

0.12 

0.13 

^^Localizaiion  of  a  hole  on  a 

two  and  four  As  neighbors  to  the  cation  site  is  represented  by  1 -center, 

2-  center  and  4-cenier  cases  respectively. 

TABLE  V.  The  calculated  lattice  di.stortion  introduecd  by  native  acceptors  in  CdGeA.ST . 


Near-neighbor  Separation^',  A 


Acceptor  center  I-center  2-center  4-center 


Hole  near  Vcd  2.38  2.50  2.56 

Hole  near  Vce  2.26  2.37  2.41 

Hole  near  Cdoe  2.97  2.71  2.62 


^  The  near-neighbor  separations  in  the  perfect  lattice  are  2.43  and  2.63  A  for  Ge-As  and  Cd-As,  respectively. 
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Atomistic  Calculations  of  Dopant  Binding  Energies 

in  ZnGeP2 

INTRODUCTION 

Zinc  germanium  phosphide  (ZnGeP2)  is  the  NLO  material  for  mid  IR  lasers  allowing 
high-power  tunability  in  the  spectral  region  of  2-5  |j.m.  An  absorption  band  around  1-2 
p.m  is,  however,  found  to  affect  the  usefulness  of  this  material.  This  band  is  attributed  to 
photoionization  of  a  deep  native  acceptor  center  associated  with  the  zinc  vacancy. 3-5 
E.xperimental  efforts  involving  selective  doping  of  the  material  are  underway  to  reduce 
the  concentration  of  this  acceptor-complex  in  the  lattice.  Earlier  work  on  the  properties  of 
dopants  in  ZnGeP2  have  been  limited  to  few  photoluminescence  and  Hall  effect  studies. 
Averkieva  et  al.  have  reported^  photoluminescence  spectra  and  hole  concentrations  for 
crystals  melt  doped  with  Cu,  Ga,  In,  Se,  Fe,  Cd  and  Si.  The  room  temperature  hole 
concentration  for  these  samples  were  reported  to  vary  from  a  low  of  2.8xlO*^/cm3  for  Si, 
nominally  the  value  .seen  in  undoped  crystals,  to  a  high  of  10 '^/cm^  for  In.  Ba.sed  on  a 
photoluminescence  study^,  it  is  reported  that  Cu  diffusion  can  increase  the  hole 
concentration  by  a  factor  of  10,000  in  ZnGeP2  by  a  proper  choice  of  the  diffusion  regime. 
Grigoreva  et  al.,  however,  have  reported*^  Hall  effect  data  on  samples  doped  with  Au,  Cu, 
Se,  Ga  and  In  and  stated  that  the  first  two  dopants  were  inactive  and  the  rest  were 
acceptors.  The  binding  energy  of  Sc  was  reported  to  be  0.40  cV  and  that  of  Ga  to  be  0.30 
eV.  In  a  more  recent  review  paper,  Rud'-^  has  indicated  that  Au,  Cu,  Ga,  In,  Se  and  Pt  are 
acceptors  with  activation  energies  in  eV  of  respectively  0.50, 0.30. 0.06  for  heavy  doping, 
0.03  for  heavy  doping  and  *0.4  for  light  doping,  0.40  and  0.50.  These  results  will  be 
compared  with  our  calculated  values  later  in  the  paper. 

We  have  initiated  an  extensive  theoretical  study5  on  ZnGeP2  and,  in  this  paper,  we  will 
report  the  results  of  atomistic  calculations  involving  cation  dopants,  namely  Cu,  Ag,  B, 
Al,  Ga  and  In.  The  choice  of  dopants  is  expected  to  indicate  the  ion-size  effect  on  the 
binding  energies  for  the  groups  I  and  III  dopants  which  lie  left  and  right  columns  of  the 
host  Zn  ion  in  the  periodic  table.  Our  approach  is  based  on  the  pair-wise  description  of 
the  lattice  which  has  been  very  successful  in  yielding  reliable  defect  energies  in  a  wide 
variety  of  materials  including  sixfold-  and  fourfold-coordinated  structures.  *0- 1 1 
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INTERATOMIC  POTENTIALS 

In  our  approach,  the  lattice  is  considered  to  consist  of  ions  interacting  via  a  combination 
of  electrostatics  and  interatomic  potentials,  while  polarization  is  included  by  means  of  the 
shell  model.*-  The  total  energy  of  a  crystal  is  then  taken  to  be  a  sum  of  two-body 
interactions  in  the  lattice.  The  pairwise  interaction  term  consists  of  the  long-range 
Coulombic  part  and  the  short-range  repulsive  part  which  is  given  by  an  analytical 
expression  of  the  Bom-Mayer  form  : 

VsR(Rij)  =  Aijexp(-Rij/pij), 

where  A  and  p  are  the  parameters  yet  to  be  determined. 

We  apply  a  parameter-free  procedure  to  model  interatomic  interactions  of  dopant  ions 
with  the  host  anion.  The  general  strategy  has  been  described  in  detail  previously**^ 
and  consist  of  two  steps  :  (i)  generation  of  quantum-mechanical  crystalline  ionic 
electronic  densities  (lEDs)  for  the  ions  involved  in  the  interaction,  and  (ii)  application  of 
the  electron  gas  formalism'-^  to  derive  the  corresponding  crystal  consistent  interatomic 
potential  (CCIP).  To  obtain  the  quantum-mechanical  de.scription  of  P-^*  embedded  in  the 
ZnGcP2  environment,  we  have  ii.sed  the  ab-initio  Perturbed  Ion  (aiPl)  model. The 
minimization  of  the  total  crystal  energy  required  by  the  Hartree-Fock  approach  provides  a 
set  of  crystal-like  atomic  wavcfunctions  that  re.spond  .self-consistcntly  to  the  nearly  exact 
crystal  potential.  In  ourca.se  the  aiPI  computation  has  been  performed  at  the  experimental 
value  of  the  lattice  parameters  of  ZnGep2  using  the  Slater  type  orbitals  given  by  Clementi 
and  Roetti.*^  It  is  interesting  to  remark  that  the  e.s.sential  electron  gas  a.ssumption  of  the 
total  crystal  density  being  a  superposition  of  the  individual  lEDs  is  better  satisfied  by  the 
aiPI  crystal-like  Orbitals  than  by  solutions  employed  in  earlier  electron  gas  calculations.' 
This  behavior  is  mainly  due  to  the  inclusion  in  the  representation  of  the  aiPI  crystal 

potential  of  a  projection  operator  that  enforces  the  ion-lattice  orthogonality. 

For  calculations  based  on  the  electron  gas  theory,  we  have  used  the  following  functional 
forms  to  evaluate  the  pairwise  interatomic  energy  r  Thomas-Fermi  for  kinetic,  Lee-Lee- 
Parr  for  exchange  and  Wigner  (with  the  new  fitting  by  Clementi)  for  correlation.  The 
coulombic  contributions  has  been  exactly  evaluated  by  means  of  analytical  algorithms 
recently  developed.  *^'*^  To  compute  the  rest  of  the  contributions,  we  have  used  a 
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.spheriodal(()),  X,  jj.)  coordinate  system.  The  angular  integration  over  (p  is  trivially  2jt.  since 
the  electron  density  is  cylindrically  symmetric  along  the  intemuclear  axis.  In  the  cases  of 
X  and  |i,  we  have  used  a  60x60  Gauss-Legendre  quadrature.  The  total  volume  of 
integration  contains  at  least  99.9998%  of  the  electron  density.  For  each  set  of  interatomic 
potential  representing  interaction  between  dopant  and  the  P  ion,  the  interatomic  energy  is 
evaluated  at  separations  starting  from  2  bohr  to  12  bohr  at  intervals  of  0. 1  bohr. 

The  calculated  potential  energy  surface  for  all  the  dopants  considered  here  is  shown  in 
Fig.  1.  The  effect  of  ion-size  is  clearly  illustrated  here;  higher  the  atomic  number  of  a 
dopant,  stronger  its  short-range  interatomic  interaction  with  P  in  the  lattice.  The  energy 
surface  is  then  fitted  to  the  analytical  expression  (Eq.  1)  whose  parameters  (i.e.  A  and  p) 
are  listed  in  Table  I.  For  the  host  lattice,  the  interatomic  potential  parameters  (i.e.  Zn-P, 
Ge-P,  P-P,  Zn-Zn,  Ge-Ge  and  Zn-Ge)  along  with  the  shell-model  parameters  Yp  and  kp 
parameters  are  obtained  by  the  empirical  fitting  using  the  program  GULP.**^  Both  dopant 
and  host  lattice  cations  are  considered  as  rigid  ions  in  the  lattice.  The  derived  interatomic 
potential  set  (Table  I)  for  the  host  lattice  reproduces  the  calculated  crystal  properties  such 
as  stnicture,  elastic  and  dielectric  constants  very  well;-*^ 

DOPANT  ENERGETICS 

The  nature  of  a  cation  dopant  in  ZnGeP2  depends  on  the  host  atom  which  it  is  replacing 
in  the  lattice.  In  the  pre.sent  case,  the  group  III  dopants  (i.e.  B.  Al,  Ga  and  In)  would  either 
act  as  a  donor  sub.stituting  Zn  or  act  as  an  acceptor  when  they  replace  Ge.  The  group  I 
cations  (i.e.  Cu  and  Ag)  would  always  be  acceptors  substituting  either  Zn  or  Ge  in  the 
lattice.  The  binding  energy  of  an  electron  or  a  hole  to  the  dopant  can  then  be  written  as  : 

AE=  Ejop+iso  ■  (Eiso  +  E(iop)  (2) 

where  Ejop+iso  refers  to  energy  of  the  dopant  complex,  Ejso  is  energy  of  an  isolated  hole 
or  electron  and  Edop  is  energy  of  the  dopant  in  the  lattice.  Calculations  to  obtain  energies 
for  Eq.  2  include,  for  example,  Cu  at  the  Zn  site  (Edop).  an  isolated  hole  localized  over 
As  ions  (Ec)  and  Cu  next  to  the  localized  hole  (Edop+c)  in  the  lattice.  Note  that  a  direct 
comparison  can  be  made  between  AE  and  the  acceptor  binding  energy  obtained  from  the 
Hall-effect  measurements. 
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Total  energies  of  defects  including  both  native  and  dopants  have  been  calculated  using 

the  Mott-Littleton  methodology-’ -22  in  which  the  lattice  is  divided  into  a  series  of 

different  regions  around  the  defect  by  concentric  spherical  boundaries.  Immediately 
surrounding  the  defect  is  region  1  in  which  all  ions  are  treated  explicitly  and  fully 
optimized.  Beyond  this  is  region  2a  in  which  all  atoms  are  still  explicitly  considered,  but 
the  relaxation  effects  are  much  smaller  and  can  be  treated  more  approximately.  The 
displacements  in  region  2a  can  be  calculated  based  on  the  force  due  to  region  1.  Beyond 
region  2a  is  region  2b  in  which  the  relaxation  energy  is  treated  yet  more  approximately 
and  responds  only  to  the  total  charge  on  the  defect  position  at  the  defect  center.  Hence 
this  component  can  be  summed  to  infinity  through  the  use  of  lattice  summation  methods 
analogous  to  the  Ewald  sum  for  the  first  inverse  power  of  distance.  In  the  present  work, 
region  1  contains  approximately  350  atoms  which  was  found  to  be  sufficient  to  converge 
the  absolute  defect  energy  to  better  than  0.01  eV. 

The  calculated  binding  energies  (AE)  for  various  dopants  are  given  in  Table  II.  At  the  Zn 
site  in  ZnGeP2,  Cu  and  Ag  acting  as  acceptors  are  predicted  to  have  the  binding  energy  of 
0.25  and  0.17  eV  respectively.  Comparison  of  the  calculated'value  with  the  experimental 
value  of  0.30  eV  for  Cu  shows  a  very  good  agreement.  The  group  III  dopants,  on  the 
other  hand,  have  a  donor-like  nature  at  the  Zn  site.  Their  electron  binding  energies  vary 
from  0.36  to  1.50  eV  with  a  significant  change  in  energy  going  from  B  to  Al.  This  may 
be  due  to  a  large  mismatch  between  the  atomic  sizes  of  B  and  Zn  in  the  lattice  (Table  II). 
For  Al,  Ga  and  In,  calculations  yield  large  binding  energies  for  the  electron  placing  the 
donor  levels  in  the  middle  of  the  gap  in  ZnGep2.  At  the  Ge  site  where  the  group  III 
dopants  act  as  acceptors,  the  results  find  a  distinct  behavior  by  B  which  forms  a  deeper 
center  as  compared  to  other  ions.  This  difference  can  again  be  understood  in  terms  of  the 
large  lattice  relaxations  arising  due  to  difference  in  covalent  radii  of  B  and  Ge  which  tend 
to  stabilize  the  bole  more  strongly.  The  group  I  dopants  (i.e.  Cu  and  Ag)  at  the  Ge  site 
have  a  very  large  hole  binding  energy  due  mainly  to  electrostatic  interactions,  for 
example,  between  [CuGe]  and  a  hole. 

In  summary,  we  have  derived  the  interatomic  potentials  for  various  cation  dopants  in 
ZnGePz  and  have  calculated  their  binding  energies  to  a  hole  or  an  electron  depending  on 
the  substitution  site  in  the  lattice.  The  calculations  based  on  atomistic  description  of  the 
lattice  predict  that  the  group-III  dopants  would  bind  the  hole  more  strongly  than  the 
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group-I  dopants.  Furthermore,  the  results  find  that  dopants  acting  as  donors  introduce 
midgap  states  in  ZnGeP2. 
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TABLE  I.  Interatomic  potential  parameters  (given  in  Eq.  2)  both  for  dopants  (obtained 
by  electron-gas  methods)  and  the  host  lattice  (obtained  by  emprical  fitting  method). 


Aij  (eV) 

CL 

electron-gas  :  Cu'+  -  P^‘ 

1562.56 

0.3694 

Agl+-  p3- 

2170.20 

0.3514 

B3+  -  p3- 

2566.54 

0.2808 

A13+  -  P3- 

2310.8 

0.3233 

Ga3+-P3- 

2567.12 

0.3474 

In3-»--p3- 

2879.71 

0.3501 

empirical :  Zn-*  -  P-3‘ 

2776.92 

0.2946 

Ge-*+  -  P-3- 

1030.93 

0.4659 

Zn2+-Zn2^- 

2776.92 

0.2946 

Ge*»+  -  Ge-»+ 

646.280 

0.1996 

Zn2+-Ge^+ 

23944.3 

0.3406 

p3-,  p3- 

1422.15 

0.4932 

4' 


% 


'j 


TABLE  II:  Binding  energies  for  various  cation  dopants  in  ZnGeP2. 


Dopant  Nature  Binding  energy*  (eV)  Rdopant**(A) 


[Cuznl 

acceptor 

0.25 

1.35 

[AgZn] 

acceptor 

0.17 

1.53 

[Bzn]' 

donor 

0.36 

0.88 

[Alzn]' 

donor 

1.05 

1.26 

(Gi^Znl 

donor 

1.40 

1.26 

(Inznl’ 

donor 

1.49- 

1.44 

[CltGel'‘‘ 

acceptor 

2.2 

1.35 

[AgCel 

acceptor 

2.0 

1.53 

lBGel\ 

acceptor 

0.83 

0.88 

[AIgcI 

acceptor 

0.61 

1.26 

[GugcI 

acceptor 

0.57 

1.26 

[ince] 

acceptor 

0.56 

1.44 

*For  Cu.  Ga  and  In  as  acceptors,  the  experimental  values^'^  are  0.30. 0.  3 1  and  0.40  eV  respectively. 
**The  covalent  radii  of  the  host  ions  ;  Rzn  :  1  3 1  A.  Rce  :  1  •22A. 


46 


* 
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First  Principles  Study  of  Stability,  Band  Structure  and  Optical 
Properties  of  the  ordered  Geo.ooSno.oO  alloy. 


Group  IV  semiconductor  alloys  have  immense  potential  for  applications  in  the  next  gen¬ 
eration  of  Si-based  electronic  and  photonic  devices.  Alloys  and  ordered  compounds  of  Si-Ge, 
Si-Sn  and  Ge-Sn  may  have  particularly  unique  optoelectronic  properties  for  applications  in 
quantum-well  intersubband  technology.^  Recently  Deng  el  al.^  have  reported  the  mor¬ 
phological  evolution  of  an  epitaxially  strained  Ge-Sn  alloy  thin  film  on  Ge  (100)  wafers.  The 
surtace  morphology  of  the  alloy  film  showed  the  formation  of  a  pattern  of  trenches  and  wires 
along  <100>  direction  created  by  the  migration  of  Sn  islands  on  the  growing  surface.  It 
is  suggested  that  the  high-throughput  fixbrication  of  novel  nanostructures  can  be  made  by 
forming  .such  trenches  on  epitaxially  strained  alloy  thin  films. 

.At  ambicmt  conditions,  both  Ge  and  fv-Sn  crystallize  in  the  diamond-like  phase  and  are 
immi.scible  in  the  l)ulk.  The  equilibrium  solid  solubility  of  Sn  in  Ge  is  le.ss  than  1%  while 
that  of  Ge  in  Sn  is  migligible.  In  spite  of  a  large  lattice  mismatch  of  about  15%  between 
Ge  and  a-tin,  epitaxially-stabilized  GeSn  thin  films  have  been  successfully  grown  on  various 
substrates.  For  example,  Hbchst  et  al:*  have  prepared  crystalline  thin  films  of  Gei_j;Sni  on 
Ge(lOO)  substrates  with  xo;0.5.  The  Sn-rich  alloy  films,  on  the  other  hand,  were  grown**  on 
InSb  and  InSb/GaAs  substrates  having  excellent  crystallinity  for  x  <  0.90  while  the  Ge-rich 
alloy  films  were  epitaxially  prepared^  on  Ge  (100)  for  x  up  to  0.24.  On  the  theoretical 
front,  the  nature  of  the  band  gap  of  these  alloys  have  been  studied®  in  the  virtual-crystal 
approximation.  The  gap  is  predicted  to  be  direct  for  0.3<x<0.75.  For  i<0.3,  the  films  are 
expected  to  exhibit  indirect  gap  while  the  films  with  x>0.75  are  expected  to  be  metallic,  as 
also  observed.® 

Calculations  based  on  first  principles  methods  are  now  routinely  performed  to  gain  an 
understanding  of  structural  and  electronic  properties  of  a  wide  variety  of  materials  including 
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the  group  I\'  (c.g.  SiGo)  binary  alloy  systems.'’'''’  Since  a  knowledge  of  structural,  electronic 
and  optical  properties  of  GeSn  alloys  is  critical  in  promoting  the  alloy  thin  films  for  device 
applications,  we  have  embarked  upon  a  theoretical  study  of  this  alloy  system  using  state-of- 
the-art  linear  combination  of  atomic  orbitals  (LCAO)  method  in  the  framework  of  density 
functional  theory  (DFT).  In  this  paper,  we  present  initial  results  on  the  ordered  Geo.soSno.so 
alloy  in  the  zincblende  phase  obtaining  its  equilibrium  structural  properties,  band  structure, 
dielectric  constant  and  reflectance  spectrum. 

In  the  LCAO  method  a  linear  combination  of  Gaussian  orbitals  are  used  to  construct 
a  localized  atomic  basis  from  which  Bloch  functions  are  constructed  by  a  further  linear 
combination  with  plane-wave  phase  factors.  The  Gaussian  basis  sets  used  here  consist  of 
five  s,  four  p  and  two  //-type  shells  for  Ge  (i.e.  a  97631/  7631/  61  set),  and  seven  s,  six  p 
auul  tliree  //-type  shells  for  Sn  (i.e.  a  9763111/  763111/  631  set).*’ 

Th(^  total  energy  cidctdations  are  performed  in  the  non-loeal  d(?nsity  approxinuition  to 
/l/msitv  fnnetional  tiuiory  (r/Tern^d  to  as  GGA)  with  a  combination  /if  Bo/.’ke  excluinge 
functional**  with  the  Perdew-VVang’s  correlation  functional.*'*  The  toh.'rance  on  the  totarl 
energy  convergence  in  the  iterative  solution  of  the  Kohn-Sham  equations  is  set  to  lO"** 
Hartrcc  and  a  grid  of  29  ^-points  was  used  in  the  irreducible  Brillouin  Zone  for  integration 
in  the  reciprocfil  space.  Under  these  computational  conditions,  the  residual  numerical  un¬ 
certainty  in  the  the.se  calctdations  is  estimated  to  be  about  0.01  eV  per  atom.  The  program 
package  CRYSTAL  is  used  for  calculations.  *‘‘  It  has  been  applied  successfully  to  study  bulk 
and  surface  properties  of  covalent  materials  such  as  Si,  C,  BN  and  GaN.*®"^^ 

We  now  proceed  with  the  geometry  optimization  of  the  GeSn  alloy  in  the  zincblende  phase 
assuming  that  structures  and  local  bonding  of  the  ordered  alloy  is  similar  to  the  constituent 
elements,  namely  Ge  and  a-Sn.  The  calculated  potential  energy  surface  of  GeSn  together 
with  its  elemental  components  Ge  and  o:-Sn  is  then  fitted  to  the  Vinet  equation  of  state 
to  obtain  the  equilibrium  structural  parameters  which  are  listed  in  Table  I.  The  calculated 
lattice  constant  follows  the  expected  trend;  the  GGA  values  are  about  1.0%  larger  than  the 
corresponding  experimental  values  for  Ge  and  a-Sn.  For  the  ordered  alloy,  calculations  find 
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4 

W  , 

A  a  (i.c.  (iceSn  -  (1/-)  +  «5n))  f<>  hc  almost  zero  obeying  Vegards  linear  rule. 

The  calculated  bulk  modulu-s  is  53  GPa  for  Geo.r>oSno..ii)  while  that  of  Ge  and  a-Sn  is  71 
and  42  GPa,  respectively  (Table  I).  The  e.xperimental  values  of  the  bulk  modulus  for  Ge  and 
ft-Sn  are  77  and  42.6  GPa,  respectively.'^  The  calculated  bulk  modulus  of  the  ordered  alloy 
therefore  appears  to  fall  on  the  straight-line  average  of  the  constituent  elements.  The  differ¬ 
ence  (AB)  between  the  alloy  bulk  modulus  and  the  average  value  Baup  (i.e.  (l/2)(BGe+Bsn)) 
comes  out  to  be  about  1.5  GPa  leading  the  ratio  of  AB/Bauj  to  be  3%.  It  is  to  be  noted 
here  that  various  theoretical  calculations  reported  the  AB/Bau^  ratio  to  be  about  0-2%  in 
SiGe.*-'’-''  The  bond-strength  in  the  ordered  alloy  is  then  predicted  to  be  stiffer  than  the 
average  of  that  iji  the  elemental  components. 

VVe  now  calculate  the  plu\se  stability  of  the  cubic  GeSn  alloy  which  is  known  to  depend 
on  the  magnituile  of  the  enthalpy  of  formation  (AH).  Following  the  notations  of  \Iartins  and 
ZiingcM-'-',  AH  is  calculated  to  be  juujtitivti  with  a  value  of  0.068  eV/atom  at  zero  pressure  ami 
temperatun!.  In  contrast,  it  is  interesting  to  not(»  hen;  that  the  LC.40-GG.\  calculations 
yield  4-0.63  and  4-0.0 1  eV'/atom  for  AH  in  the  cubic  alloys  ofSn  with  C  and  Si,  respectively.'^' 
The  calculated  results  therefore  indicate;  the  relative  stability  of  the  Gebn  alloy  against 
.segregation  into  the  elemental  components.  Furthermore,  calculations  of  Gibbs  free  energy 
show  the  stability  of  the  cubic  phtxse  of  GeSn  at  presstires  up  to  9  GPa  at  which  one  of  the 
alloy  components,  Ge  undergoes  a  high-pressure  phase  transition.”- 

VVe  now  turn  our^ttention  to  electronic  properties  of  GeSn  and  display  the  upper  valence 
and  lower  conduction  band  structure  of  the  ordered  alloy  in  Fig.  1  .  Accordingly,  the  top 
of  the  valence  band  consists  of  the  triply  degenerate  Fis  in  the  absence  of  the  spin-orbit 
interaction  terms  in  this  work.  The  p  states  of  Ge  and  Sn  forms  the  upper  valence  band 
which  has  a  width  of  about  6.3  eV.  The  s  states  of  Ge  and  Sn  form  a  band  with  a  band-width 
of  1.6  eV  at  about  10  eV  below  the  top  of  the  valence  band.  The  Sn  4d  band  appears  at 
about  23  eV  below  the  top  of  the  valence  band  and  is  therefore  least  likely  to  influence  the 
states  in  the  conduction  band  (Fig.  2).  In  the  ordered  alloy,  dispersion  of  the  valence  band  is 
similar  to  that  of  Ge  in  agreement  with  the  photoemission  study.®  However,  its  conduction 
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band  appears  to  bo  similar  to  that  of  a-Sn  in  which  tin?  minimum  occurs  at  f.  The  valence 
band  offsets  in  the  alloy  components  are  calculated  to  be  -0.31eV  (Ge)  and  -0.r2eV  (Sn) 
with  respect  to  the  reference  energy  of  the  top  of  the  valence  band  of  the  cubic  GeSn. 

The  lowest  gap  between  top  of  the  valence  band  to  bottom  of  the  conduction  band  is 
found  to  be  at  F  in  the  ordered  alloy.  Thus,  the  present  DFT  calculations  predict  the  gap  to 
be  direct  with  a  gap  value  of  0.1  eV.  Knowing  that  the  DFT  calculations  underestimate  the 
band  gap  relative  to  experiments,  we  now  estimate  the  gap  from  the  experimental  gap  values 
of  the  alloy  components.  cv-Sn  is  a  aero-band  gap  semiconductor  while  Ge  is  an  indirect-gap 
semiconductor.*^’"'*  The  minimum  indirect  gap  in  Ge  is  0.78  eV  associated  with  the  F-L 
transition  while  the  minimum  direct  gap  is  0.98  eV  at  F.  Taking  the  theoretical  value  ot 
the  optical  bowing  parameter  of  0.8eV  at  F,  the  direct  gap  (at  F)  in  the  alloy  is  estimated*** 
to  be  0.3  eV.  It  is  to  be  noted  here  that  angle-resolved  photoemission  study  on  Geo..i8Sno.52 
thin  films  does  find  a  narrowing'of  the  alloy  band  gap  with  respect  to  the  average  of  the  gap 
of  its  components.'* 

Both  valence  and  conductioji  level  energies  at  F,  -V  and  L  k-pointsshow  a  linear  variation 
with  applied  pressure.  The  calculated  pressure  coefficients  of  the.se  levels  listed  in  Table  II 
yield  the  coefficient  of  the  minimum-energy  direct  gap  at  F  to  be  3.45  x  lO"**  cV/bar.  At 
about  7  GPa,  a  cro.s.s-over  between  the  conduction  levels  at  F  and  A'  occurs  resulting  into 
the  minimum-energy  gap  to  be  indirect  in  GeSn. 

It  is  well  known  that  the  conduction  bands  are  not  well  described  by  the  density  functional 
theory,  though  their  use  in  polarizability  calculations  yield  satisfactory  values  of  dielectric 
response  functions  in  semiconductors.*®'^**  In  this  work,  a  siim  over  states  (SOS)  method 
is  used  to  calculate  the  real  and  imaginary  parts  of  the  polarizability  from  which  related 
functions  such  as  the  dielectric  constant  and  energy  loss  function  (ELF)  can  easily  be  ob¬ 
tained.  The  SOS  calculations*®-^®  yield  the  value  of  44.1  A®  for  (static)  polarizability,  19.6 
for  dielectric  constant  and  4.36  for  index  of  refraction  of  the  ordered  alloy.  At  6889  nm, 
the  index  of  refraction  is  calculated  to  be  4.44  which  is  very  close  to  the  Soref’s  empirical 
estimation*  of  4.52.  It  is  to  be  noted  here  that  the  calculated  dielectric  constant  of  16.0  and 
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2-1.3  for  G(?  iuul  a-Sn  arc  arc  in  very  good  agreement  with  the  reported  values"'*  of  16.0  and 


23.8. 

The  complex  dielectric  response  function  is  now  used  to  calculate  the  plasmon  energy 
corresponding  to  the  maximum  ot  the  ELF  hinction.  It  is  found  to  be  16.1  e\  for  GeSn 
and  is  smaller  than  the  average  of  the  component  values  of  19.3  eV  and  14.7  eV  for  Ge 
and  a-Sn,  respectively.  Our  calculated  plasmon  energy  for  Ge  can  be  compared  with  the 


previously  reported  value*'^  of  about  1  <  .5  eV  where  the  real  part  of  the  dielectric  constant 
was  obtained  from  the  imaginary  one  via  the  Kramers-Kronig  relation.  The  reflectance, 
stiuare  of  the  normal  incidence  reflectivity  is  determined  from  the  dielectric  function  via  the 


Fresnel  equations,  with  respect  to  the  incident  photon  energy  for  GeSn  along  with  Ge  and 
Sn.  As  expected,  the  reflectance  behaviour  of  the  ordered  alloy  can  be  approximated  by  the 
average  of  the  reflectance  of  Ge  and  Sn.  Similarly,  the  mean  value  of  the  Compton  profile^^ 


(with  respect  to  the  moment  v;du(^  (j  of  electron  .in  the  scattering  direction)  ot  GeSn  i 


'A  rv: 


to  l)(^  the  average  of  that  in  the  ehunental  components. 

In  summary,  first  principles  calculations  basc'd  on  tlu'  density  functional  theory  find  the 
ordered  GeSn  alloy  to  be  stal)l(‘  in  tlu?  zincbh^iuh?  ph.use  up  to  9  GPa.  The  band  gap  of 


tlu;  cubic  allov  is  pr(;dicte<l  to  lu;  din;ct.  Although  the  bond-strength  in  the  ordered  allov 
is  found  to  be  stiffer  than  the  average  ot  that  in  the  elemental  components,  the  reflectance 
and  Compton  profile  of  the  alloy  can  be  obtained  by  taking  the  average  of  their  respective 
values  in  the  elemental  components. 
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sum  of  oscillator  strengths  to  be  equal  to  the  number  of  electrons  per  unit  cell. 

C.  S.  Wang  and  B.  M.  Klein,  P/ij/s.  Reu.  B  24,  3417  (1981). 
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TABLES 


TABLE  I.  Structural  properties  of  Geo.soSuo.sOi  Ge  anti  cv-Su. 


a(A) 

V(A3) 

B(GPa) 

B' 

GeSn 

This  work  (GGA) 

6.20 

59.58 

53 

6.9 

Ge 

This  work  (GGA) 

5.74 

47.28 

71 

5.1 

Expt.(Ref.  24) 

5.66 

45.33 

77 

— 

tK-Sn 

This  work  (GGA) 

6.65 

73.38 

42 

5.1 

Exi)t.(Rcf.  19) 

6.49 

68..34 

•  42.6  . 

-- 

TABLE  II. 

Electronic  pro|)crti«\s  of  Gt^o.noSno.rjo  in  tin?  '/inchltnulo  plnuso. 

Danr!  Gap,  eV 

Direct 

r 

O.l 

X 

3.1 

L 

1.6 

Indirect 

r-  X 

1.0 

r  -  L 

0.7 

Pressure  Coefficients,  x  10“** 

eV/bar 

V^alcnce  Band 

Tis 

2.40 

Xs 

1.13 

L3 

1.85 

Conduction  Band 

r. 

5.85 

Xi 

1.70 

Li 

3.77 

00 


